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FOREWORD 


Hiis  final  report  was  prepared  by  Aeronautical  Research 
Associates  of  Princeton,  Inc.,  Princeton,  New  Jersey,  in  fulfill¬ 
ment  of  USAF  Contract  Number  AF33(6l5)-2l44;  mm%  ^6899-6x430014- 
0000-600-FD);  Project  and  Task  Numbers:  None.  Bits  effort  was  done 
under  Laboratory  Director's  Funds.  This  project  is  part  of  a 
program  to  provide  the  theoretical  and  practical  applied  research 
necessary  to  maintain  USAF  structural  design  criteria  for  aerospace 
flight  vehicles  and  to  establish  optimum  structural  requirements, 
flight  limitations  and  design  philosophies  commensurate  with 
vehicle  structural  life. 

The  work  was  administered  under  the  direction  of  the  Air  Force 
Flight  Dynamics  Laboratory,  Research  and  Technology  Division,  Air 
Force  Systems  Command,  Weight -Patterson  Air  Force  Base,  Ohio,  Mr. 
Paul  L.  Hasty  (FDTR),  Project  Engineer.  This  investigation  was 
administered  under  Project  1367  and  Task  136702;  all  follow-on  work 
will  be  funded  under  this  task. 

The  work  reported  in  this  study  was  conducted  by  Aeronautical 
Research  Associates  of  Princeton,  Inc.  with  Dr.  John  C.  Houbolt  as 
principal  investigator,  and  covers  the  period  from  November  1964  to 
March  1966.  The  manuscript  was  released  by  the  author  in  March  1966 
for  publication  as  an  RTD  Technical  Report. 

The  Contractor's  Report  Number  is  ARAP  Report  No.  83,  Volumes 
I  and  II. 


ABSTRACT 


A  review  is  made  of  the  progress  and  findings  that  have  been 
obtained  in  a  study  aimed  at  the  development  of  gust  design  proce¬ 
dures  based  on  power  spectral  techniques.  Included  is  a  consider¬ 
ation  of  atmospheric  turbulence  makeup.  A  corresponding  discussion 
is  given  of  the  aircraft  structural  parameters  that  are  needed  in 
spectral  evaluation  and  on  the  practical  aspects  of  their  evalua¬ 
tion.  Emphasis  is  given  to  the  spectral  design  procedures  that 
have  evolved  from  the  study,  and  a  new  hypothesis  for  fatigue 
considerations  is  advanced,  cast  in  the  notation  of  gust  loads 
design.  Ease  of  application.  Interpretation  and  implication  to 
design,  present  and  future,  and  comparison  with  past  design  proce¬ 
dures  are  covered. 
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SYMBOLS 


&  lift  curve  slope 

A  structural  parameter;  <?x  *  Ao w 

c,cQ  wing  chord 

f  frequency,  cps 

h  altitude 

H(cd)  frequency  response  function  due  to  sinusoidal  gust 

k  reduced  frequency,  k  *  oe/2V 

K  gust  alleviation  factor 

g 

L  scale  of  turbulence 

n  total  number  of  times  load  level  x  Is  crossed  with 

positive  slope 

nQ  total  number  of  zero  crossings  with  positive  slope 

N  number  of  times  per  second  load  level  x  Is  crossed  with 

positive  slope 

Nq  number  of  zero  crossings  per  second  with  positive  slope 
p(o)  probability  distribution  function  of  a 

P  proportion  of  time  spent  in  turbulence 

sr  reference  stress  as  used  in  fatigue 

S  wing  area 

T  total  (lifetime)  flight  time 

Tt  total  time  in  turbulence 

Tp  average  duration  of  a  mission  flight 
Ud  gust  velocity  value  similar  to  Ude 
Ude  "derived"  gust  velocity 

V  flight  velocity 

Ve  equivalent  airspeed 


SYMBOLS  (Cont'd6) 


W  airplane  weight 

x  Increment  In  load  or  response  level  due  to  gusts 

Tj  load  level  representing  limit  load 

x,  1  -  g  load  level 

a, (3  constants 

X  wavelength 

u  mass  ratio  or  mass  parameter,  2W/apegS 

o 

M-q  mass  ratio,  |1q  -  %g 

p  air  density 

a  r.m.s.  value 

cr  composite  r.m.s.  value  of  vertical  gust  velocity 
a  r.m.s.  value  of  response  x 

A 

a  r.m.s.  value  of  vertical  gust  velocity 

<t>(Cl)  power  spectrum 

cd  angular  frequency 

Cl  spacial  frequency,  cd/V 

Note: 

Cl  =  *  =  ^  =  2L 

VOX 

4>(Cl)  =  =  v<£>(cd)  =  |  0(k) 
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SECTION  1 


INTRODUCTION 

This  report  deals  with  the  subject  of  designing  aircraft  for 
gust  loads  encounter  on  a  power  spectral  basis.  Recent  studies 
suggest  that  developments  on  the  subject  during  the  past  several 
years  could  be  simplified  or  streamlined  considerably  and,  in  fact, 
that  procedures  could  be  put  on  a  much  firmer  basis,  while  still 
preserving  generality.  It  is  the  purpose  of  this  report  to  bring 
out  for  consideration  the  results  of  these  newer  studies,  and  to 
indicate  the  associated  improvements  that  are  desired  in  the 
analysis  of  flight  operational  data.  At  the  same  time,  a  new 
hypothesis  is  advanced  for  the  treatment  of  fatigue  due  to  gust 
loading  that  circumvents  the  necessity  of  introducing  various 
notions  such  as  cumulative  damage  concepts  and  the  like.  The 
procedure  for  treating  the  fatigue  problem  simultaneously  -with  the 
gust  loads  exceedance  problem  is  shown.,  uBing  a  common  language 
and  notation  for  both. 

The  body  of  the  report  presents  the  essential  results,  while 
appendices  are  used  for  the  more  detailed  treatments.  Appendix  A 
defines  and  indicates  some  of  the  practical  aspects  that  are 
involved  in  the  determination  of  the  basic  structural  response 
parameters  A  and  NQ  ;  Appendix  B  gives  a  brief  development  of 
the  gust  design  treatment  advanced  here,  and  Appendix  C  covers  the 
new  fatigue  hypothesis. 
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SECTION  2 

BASIC  RELATIONS  FOR  LOAD  COUNT 

Oust  Encounter 

As  a  result  of  examining  past  theoretical  studies  and  airplane 
operation  statistics  on  gu3t  loads,  and  on  the  basis  of  further 
theoretical  study  which  is  reviewed  briefly  in  Appendix  B,  a  newer 
and  concise  concept  on  the  statistical  treatment  of  gust  loads  has 
been  developed.  The  concept  is  embodied  in  a  single  universal 
curve.  Figure  1,  which  allows  for  the  rather  quick  estimation  of 
the  expected  number  N  of  times  per  second  given  load  levels  are 
crossed  during  flight  of  an  airplane.  The  definition  and  meaning 
of  the  other  terms  follow. 

A  and  Nos  These  are  the  basic  airplane  structural  response 
parameters  which  apply;  A  relates  in  linear  fashion  the  r.m.s. 
value  a  of  gust  ?nput  to  the  r.m.s.  value  a  for  the  response 

v  A 

according  to  the  relatior 

=  Aoc 

The  parameter  Nq  is  a  characteristic  frequency  parameter  which 
denotes  the  number  of  times  per  second  the  response  time  history 
crosses  ^e  value  zero  with  positive  3lope.  Appendices  A  and  B 
define  these  parameters  more  completely  and  indicate  the  practical 
aspects  of  their  evaluation. 

x:  x  refers  to  the  chosen  response  quantity,  sucn  as  bend¬ 
ing  stress,  bending  moment,  acceleration,  etc.  More  precisely, 
it  is  an  incremental  value  Ax  above  the  1  -  g  load  level  value, 
and  for  design  is  expressed  in  terms  of  chosen  design  allowables 
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in  accordance  with  the  design  philosophy  being  followed;  fo: 
example,  for  design  in  terms  of  limit  load,  x  is 


x  =  Ax  =  xL  L 


where  x^  ^  is  the  allowed  limit  load  value. 

Qq  and  P:  a  is  a  composite  value  denoting  the  r.m.s. 
value  of  gust  velocity  for  all  turbulence  encounter  during  the 
lifetime  of  the  aircraft,  and  P  is  the  associated  value  denoting 
the  proportion  of  total  flight  time  that  the  airplane  is  in  turbu¬ 
lence;  these  will  be  discussed  in  greater  detail  in  the  section  on 
design. 

The  curve  of  Figure  1  is  found  (Appendix  B)  to  be  represented 
well  by 


_N_  /3.36  Asv6-?4 

PN0  \  x  ) 


(1) 


for  x/A a  >  9.13;  a  relation  applying  to  lower  values  of  the 

V 

abscissa  is  given  in  Appendix  B  as  well,  but  since  this  range  is 
not  normally  of  interest  the  relation  is  not  repeated  here.  The 
empirical  constants  3.36  and  6.74  are  tentative  in  nature,  and 
should  be  updated  as  more  and  more  operational  data  are  cast  in 
the  "universal"  form  suggested  here.  The  total  number  n  of 
times  given  load  levels  are  crossed  or  exceeded  during  the  over-all 
flight  time  of  the  aircraft  follows  simply  from  Equation  (l)  as 

6.74 


/3  .36  Aa  n 
PTNo(— T^J 


(2) 


where  T  represents  the  total  flight  time.  As  shown  subsequently, 
this  relation  forms  the  basis  for  the  design  for  gusts  from  a 
number  of  load  exceedances  point  of  view. 

Fatigue 

At  the  same  time,  because  of  the  implications  of  recent  fatigue 
data  for  random  loading  conditions,  and  because  the  manner  of  hand- 
ling  and  presenting  data  for  both  fatigue  and  gust  load  count  can 
be  developed  along  closely  similar  lines,  a  new  fatigue  concept  is 
advanced  here.  The  hypothesis  is:  fatigue  life  is  a  function 
simply  of  the  r.m.s.  stress  level  of  the  over-all  time  history  of 
stress  exposure,  regardless  of  the  sequence  of  load  application. 

In  terms  of  the  gust  loads  parameters,  the  treatment  on  fatigue 
given  in  Appendix  C  indicates  that  fatigue  life  may  be  expressed 
by  an  equation  of  the  form 


where  ng  refers  to  the  total  number  of  times  the  zero  stress 
level  may  be  crossed  with  positive  slope  before  fatigue  failure 
(on  the  basis  of  zero  mean  value),  and  where  a  and  (3  are 
empirical  factors  which  are  associated  with  the  material  and  type 
of  construction  practices  used.  The  product  £x  is  a  reference 
empirical  stress  of  the  type  often  found  in  analytical  representa¬ 
tions  of  fatigue  life,  and  hence  £  relates  this  reference  stress 
to  the  chosen  design  value  of  x;  the  product  Ac  is  noted  to  be 
simply  the  r.m.s.  value  a  of  the  random  stress  history. 

A 

The  remarkable  correspondence  between  Equations  (2)  and  (3) 
will  be  dealt  with  in  the  next  section. 
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SECTION  3 
DESIGN 


0c  and  P  Values 

The  application  of  Equation  (2)  requires  that  consideration 
be  given  to  the  choice  of  a  and  P  values.  Two  procedures 
are  offered.  In  one,  a  and  P  are  simply  stipulated  constant 

values,  and  hence  the  procedure  is  very  much  akin  to  the  discrete 
gust  approach  in  which  design  gust  velocity  Udg  of  50  fps,  for 
example,  is  stipulated  and  no  account  is  taken  of  the  amount  of 
turbulence  encounter.  This  first  procedure  is  felt  to  have 
considerable  merit  due  to  its  simplicity. 

/ 

In  the  second  procedure,  explicit  consideration  is  given  to 
the  various  turbulence  intensity  levels  encountered  and  to  the 
amount  of  time  spent  in  turbulence  at  various  altitudes.  Many 
studies  along  this  line  have  been  persued  in  recent  years.  They 
generally  have  proceeded  by  breaking  up  the  turbulence  encounter 
into  severity  brackets,  usually  of  two  types,  nonstorm  and  storm. 
Operational  data  are  used  in  an  attempt  to  deduce  the  statistically 
applicable  variations  with  altitude  of  intensity  and  proportion  of 
time  for  each  of  the  severity  brackets  chosen,  with  the  object  of 
using  these  deduced  values  later  for  design  purposes.  It  is  felt 
that  bi’eakdowns  of  this  type  are  perhaps  an  over  refinement  or  that 
the  accuracy  of  the  results  available  does  not  warrant  this  much 
detail.  Here  the  use  of  single  measures  for  gust  intensity  and 
proportion  of  time  is  advocated;  specifically,  in  question  form, 
what  statistically  is  the  over-all  r.m.s.  gust  velocity  value  of 
the  turbulence  encountered  at  each  altitude,  and  likewise,  what  is 
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the  proportion  of  time  on  the  average  that  turbulence  is  encountered 
at  each  altitude? 

Figure  2  shows  some  tentative  initial  variations  of  cl  and 

c 

P  that  are  estimated  from  the  consideration  of  existing  data  and 
analyses.  Most  of  the  data  indicate  that  ac  is  nearly  constant 
with  altitude  and  in  the  neighborhood  of  3  to  3.5  fps;  a  constant 
value  of  3.0  has  been  chosen  here.  For  P  the  choice  is  more  in 
doubt;  some  results  suggest  curve  a,  others  curve  b.  (The  circum¬ 
stance  here  is  quite  odd,  since  of  all  quantities  P  should  be 
the  easiest  to  establish.)  It  is  advocated  therefore  that  the 
appropriate  P  variation  curve  (as  well  as  the  a  curve)  be 

V 

established  from  the  data  that  exist  and  from  those  being  gathered. 
At  the  moment,  curve  a  is  felt  to  be  more  appropriate,  although 
proof  of  thi3  feeling  is  lacking. 

Thus  even  in  the  more  involved  second  procedure  discussed 
here,  only  the  P  variation  has  to  be  considered,  since  a 

c 

appears  constant.  The  effective  value  for  P  applying  to  chosen 
missions,  or  to  flight  paths  which  are  followed  on  the  average, 
can  quite  easily  be  found  from  the  P  variation  curve,  by  appro¬ 
priately  weighing  the  times  at  each  altitude  bracket.  Figure  3 
shows  example  results  for  the  assumed  flight  plan  variation  shown 
also  in  the  figure,  assuming  curve  b  of  Figure  2  is  the  appropri¬ 
ate  P  variation;  hold  times  can  be  taken  into  account  convenient¬ 
ly  by  using  an  effective  rate  of  climb  or  descent,  and  the  rather 
low  value  of  300  feet  per  minute  used  here  reflects  this  consider¬ 
ation.  Hence,  if  an  airplane  is  to  be  designed  to  fulfill  an 
average  mission  covered  by  this  figure,  then  the  effective  P  to 
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be  used  in  design  can  be  ^aken  directly  from  the  figure. 

Direct  Design 

The  direct  design  for  gust  loads  on  a  load  count  basis  and 
for  fatigue  considerations  as  well,  follow  in  a  straightforward 
way  from  Equations  (2)  and  (3).  These  equations  plot  as  shown  in 
Figure  4.  The  use  of  Aa^x  versus  PTNq  is  chosen  simply  on 
the  basis  of  the  close  resemblance  to  the  familiar  S-N  diagrams  of 
fatigue.  Suppose  that  x  is  associated  with  limit  load  values 
and  that  n  is  stipulated  not  to  exceed  5»  A  design  is  checked 
then  by  entering  the  chart  at  the  computed  values  of  Aa^x  and 
PTNq,  and  if  the  point  falls  below  the  n  =  5  line,  the  design 
is  judged  safe.  (Note,  consideration  of  the  curve  of  Figure  1  at 
the  larger  abscissa  value  indicates  that,  for  x  on  a  limit  load 
basis,  an  n  in  the  neighborhood  of  50  corresponds  to  about  one 
encounter  for  x  presented  in  terms  of  ultimate  load.) 

At  the  same  time,  safety  from  fatigue  failure  can  be  checked. 
The  particular  material  and  construction  used  fix  a  and  0, 
and  hence  fix  the  fatigue  curve.  The  remoteness  of  the  design 
point  indicates  then  the  margin  of  safety  against  fatigue  failure. 
Note  that  although  not  explicitly  stated  here,  the  fatigue  curves 
used  should  reflect  the  effect  of  the  1-g  mean  flight  stresses. 

Thus,  this  single  and  concise  chart  embodies  the  complete 
procedure  advanced  herein  for  designing  aircraft  for  gust  encounter, 
and  allows  the  nearness  to  fatigue  troubles  due  to  gusts  to  be 


monitored  as  well. 


Load 


ati 


An  illustrative  example  is  given  here  to  indicate  the 
approach: 

Design  conditions  given 

x  =  28,000  psl 
13  =  2.5 

T  =  30,000  hrs  =  10.8  X  107  sec 

n  =  5 

For  the  mission 
Oc  =  3  fps 

Tp  =  3  hrs  1  p  =  .085  from  Figure  3 
h  =  20,000  ftj 

From  the  airplane  preliminary  design 

A  =  310  psi/fps 
Nq  =  1.1  cps 

These  values  yield 

Aac 

IT  ■  -0332 

PTN0  =  10.1  X  106 

which  spots  Just  below  the  n  =  5  curve  of  Figure  4,  thus  indicat¬ 
ing  a  safe  design.  The  point  is  also  noted  to  be  far  removed 
from  the  fatigue  curves . 

Comparison  Design 

An  even  simpler  design  procedure  may  be  used  by  a  comparison 
technique.  Write  Equation  (2)  in  terms  of  airplane  2,  and  also  in 
terms  of  airplane  1,  and  divide  the  equations.  The  result  is  that 
shown  in  Figure  3.  Suppose  that  airplane  1  is  an  older  gust-eriti- 
cal  airplane  that  has  proven  itself  airworthy  by  many  successful 
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Figure  5.  Gust  Design  by  Comparison  (Note,  stresswise 
adding  metal  lowers  A). 


a  Increase  Due  to  Extended  Severe  Turbulence  Encounter 
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years  of  flight,  and  suppose  airplane  2  is  the  new  airplane  under 
design  consideration.  The  ratios  of  the  Aoc/x  and  PTNq  aua$>£3  - 
ties  ar$  established  and  the  adequacy  of  the  new  design  Is  sSmpl” 
established  according  to  whether  the  point  falls  on  the  safe  side 
or  the  unsafe  side  of  the  line.  For  airplanes  of  similar  category, 
this  appears  to  be  an  attractive  approach  and  has  the  merit  of 
good  simplicity.  Other  advantages  are  also  offered;  for  example, 
the  use  of  ratios  tends  to  minimize  uncertainties  that  may  exist 
in  the  determination  o^  the  parameters  such  as  A  and  Nq,  in 
that  errors  common  to  each  determination,  such  as,  say,  in  the 
choice  of  lift  curve  slope,  tend  to  cancel  one  another. 
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SECTION  4 

CONCLUSIONS  AND  RECOMMENDATIONS 


Two  points  are  discussed  here.  One  is  in  reference  to  the 
an  and  F  values  that  are  chosen  for  design.  The  values  and 
curves  presented  apply  basically  to  medium  and  higher  altitude 
flights  and  for  conditions  on  the  average  over  long  periods  of 
flight  time.  For  special  purpose  missions,  such  as  low  altitude 
flight  in  a  high  percentage  of  severe  type  turbulence,  the  values 
of  a  and  P  must  be  selected  accordingly.  Figure  6  is  illus- 

V 

trative  of  how  might  increase  when  operating  under  more  severe 

turbulence  conditions!  the  majority  of  the  flight  is  considered  to 
occur  at  the  r.m.s.  level  and  effectively  a  certain  percent- 

age  of  the  more  severe  level  a 2  is  added;  +  ?2  becomes  the 
total  P.  As  an  example.  If  is  taken  as  3,  then  for  the 

specific  case  of  P^/P i  =  ,05  and  og  =  7.5  the  effective  a, 
is  found  to  be  3-72. 

The  second  point  is  relative  to  the  scale  of  turbulence  L 
which  appears  in  the  expression  representing  the  power  spectrum 
for  the  input  gust  velocities.  Equation  (Al)  in  Appendix  A.  This 
quantity  was  not  touched  upon  In  the  preceding  discussion,  but 
comment  3s  merited.  In  general,  the  structural  parameter  A  13 
found  to  be  dependent  on  L,  and  therefore  the  value  of  L 
selected  for  use  is  quite  important.  The  evaluation  of  L  has 
been  of  much  concern,  and  estimates  range  all  the  way  from  500  to 
5,000  feet.  The  primary  reason  why  its  evaluation  is  illusive  at 
the  moment  is  due  to  the  fact  that  the  time  history  records  of 
vertical  gust  velocity  that  have  been  deduced  from  flight  studies 
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generally  show  uncertain  or  unreliable  low  frequency — or  long 
wavelength — content.  The  long  wavelength  components  present 
lead  to  a  large  evaluated  scale  value.  It  is  suspected,  however, 
that  these  longer  wavelengths  arise  because  of  instrumentation 
limitation  or  inaccuracies,  and  thus  the  scale  values  deduced  are 
in  question.  The  essential  point  to  be  made  here  is  that  more 
attention  should  be  paid  to  the  instrumentation  used  in  flight 
tests  and  what  can  be  expected  of  it.  A  simple  evaluation  proce¬ 
dure  that  should  be  used,  for  example,  but  which  does  not  seem  to 
have  been  followed,  is  the  following.  Plight  tests  of  the  instru¬ 
mented  airplane  should  be  made  in  still  air,  both  for  the  condi¬ 
tions  of  level  flight  and  with  deliberate  pull-up  and  push-down 
maneuvers;  these  flights  should  be  of  at  least  4  minutes  duration. 
The  records  should  then  be  evaluated  in  the  same  manner  as  if 
turbulence  had  been  encountered,  to  see  if  the  vertical  velocity 
time  history  evaluates  to  zero,  as  it  should.  Any  residual  trace 
that  appears  serves  to  indicate  instrumentation  limitations,  and 
possibly  what  low  frequency  contaminations  are  likely  to  be  found 
in  the  gust  records.  A  guide  would  then  be  available  to  select 
appropriate  low  frequency  filtering  devices  in  the  record  evalua¬ 
tion  procedures. 
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APPENDIX  A 


STRUCTURAL  PARAMETERS 


This  appendix  serves  simply  to  review  the  definition  of  two 
basic  structural  parameters,  A  and  Nq,  that  are  significant 
in  the  spectral  determination  of  gust  loads  (see  Reference  1  for 
a  fuller  treatment).  The  basic  relation  between  the  input  gust 
spectrum  and  the  output  response  spectrum  is 


0X( cd)  =  |h(cd)  | ^  0w(cd) 


where  x  is  any  dynamic  response  quantity  of  concern  and  H  is 

the  frequency  response  function  for  x  due  to  unit  sinusoidal 

gusts.  For  0  ,  an  expression  of  growing  favor  and  which  exhibits 
w 


an  (D 


-5/3 


fall  off  at  large  frequencies  is 

o^L  1  +  |(1.339  W)2 

*„(«)  -  — - 1 


W 


r  tU 

[i  +  (1.339  Ln)2J6 


where  a ,  the  r.m.s.  gust  velocity,  and  L,  the  turbulence  scale, 
characterize  the  turbulence  being  encountered,  and  where  =  a>/V. 

The  basic  structural  parameters  are  derived  from  the  output 
spectrum  in  accordance  with  the  following  sketch  and  equations. 
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The  parameter  A  relates  the  r.m.s.  value  of  gust  input  in  linear 
fashion  to  the  r.m.s.  value  of  output  according  to  the  relation 

ax  =  Aaw  (Al) 

and  Nq  Is  an  over-all  characteristic  frequency  of  response  of 
the  airplane  and  specifically  denotes  the  number  of  time3  per 
second  that  the  response  time  history  crosses  the  value  zero  with 
positive  slope. 

Instead  of  using  a  theoretical  upper  limit  of  infinity  in  the 
integration,  a  cut  off  frequency  co  is  Introduced  in  practice, 

V 

otherwise  Nq  would  appear  unrealistically  unbounded  in  some 
cases.  The  value  of  <£>c  to  use  in  practice  is  of  much  concern, 
but  a  ale  of  thumb  that  holds  much  promise  is  that  brought  out 
in  the  following  sketches  showing  the  variation  of  A  and  Nq 

with  03  . 
c 
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In  the  determination  of  the  parameter  A  there  is  not  much 
difficulty,  for  a  plateau  is  usually  reached  quite  sharply.  For 
NQ,  no  difficulty  Is  encountered  in  some  cases  such  as  case  a 
Illustrated.  For  case  b.  however,  the  difficulty  of  ever  Increas¬ 
ing  Nq  arises;  the  rule  of  thumb  here  Is  to  take  as  a>c  the 
point  where  A  has  become  essentially  flat,  point  1,  and  read 
off  Nq  at  this  frequency,  point  2.  Thus,  since  the  contribution 
to  the  total  power  or  to  A  of  spectral  information  beyond  co 

V 

Is  essentially  negligible,  the  contribution  to  Nq  beyond  o>c  is 
essentially  hash  or  noise  type  Information,  and  hence  should  be 
disregarded  In  arriving  at  level  crossing  which  has  meaning  from 
a  practical  point  of  view. 


APPENDIX  B 


LOAD  LEVELS  DUE  TO  GUST  ENCOUNTER 

We  derive  here  a  fundamental  curve  for  load  level  crossings 
that  forms  the  basis  for  the  complete  design  of  aircraft  due  to 
gusts  on  a  power  spectral  basis.  A  remarkable  closeness  to  fatigue 
strength  considerations  appears  and  so  it  is  fitting  to  advance  a 
new  and  analogous  procedure  for  fatigue  design  as  well. 

We  start  from  a  slightly  different  viewpoint  from  that  con¬ 
sidered  in  other  studies,  such  as  in  Reference  1,  but  show  that 
similar  end  results  are  obtained.  Consider  the  timewise  variation 
of  gust  encounter  of  an  airplane  during  its  lifetime.  A  time 
history  record  of  the  gust  velocities  encountered,  or  of  some 
resulting  structural  response  quantity,  might  appear  for  example 
as  shown  in  the  following  sketch 


— — 


"'Wr11 


For  convenience,  let  us  now  lump  in  link  fashion  the  quiescent 
periods  together,  and  likewise  the  turbulence  patches,  so  that  an 
equivalent  time  history  record  might  appear 


T 
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As  has  been  verified  in  a  gross  sense,  let  us  assume  that  each 
patch  of  turbulence  encountered  is  locally  stationary  and  Gaussian 
in  character.  Then  from  the  relation  due  to  Rice  (Reference  2) 


N  =  NQe 


giving  the  number  of  times  per  second  that  a  given  level  x  is 

O 

crossed  with  positive  slope,  and  with  Equation  (Al),  we  can  write 
the  number  of  time3  per  second  a  given  load  or  response  level  is 
crossed  as  the  integrated  result 


(Bl) 

r 


T. 


N 


■  kf  N°e 


o*2  2 

2A  a 


dt 


where  Nq,  A,  and  o  are  in  general  functions  of  time.  Note 
a  refers  here  to  the  r.m.s.  value  of  gust  velocity  a  .  but  for 
brevity  in  notation  we  drop  the  subscript  w.  Now  Nq  varies 
only  slightly  with  time  and  so  we  regard  it  constant.  For  general 
ity,  we  could  consider  the  timewise  variation  of  the  product  A o, 
but  since  the  variation  of  A  is  not  large,  we  also  choose  A 
constant;  specifically,  the  r.m.s.  value  A  is  suggested  as  the 

V 

logical  choice  in  accordance  with  the  approximation 

J  A2o2dt  =  A^  Jo2 dt 
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Equation  (Bl)  thus  may  be  written 


2  2 
2Ac°c 


(B2) 


where  P  denotes  the  proportion  of  total  flight  time  spent  in 
discernible  turbulence  and  oc  is  the  r.m.s.  value  of  a, 
defined  by 


The  oc  value  may  also  be  written 


=  J  a2p(a) 
0 


da 


where  p(a)  represents  the  probability  distribution  function  of 
the  r.m.s.  value  c  of  gust  velocity  that  applies  during  the 
entire  flight  experience.  These  equivalent  relations  imply  then 
that 

d(~-)  =  p(a)da  (B3) 

It 
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or 


£-/*’< °>d0 - 1 
r  0 


and  that  for  the  limits 


t  =  0, 

t  =  Tt, 


a  «  0 

O  =  oo 


With  Equation  (33),  Equation  (B2)  may  be  expressed  as 


2 

x  c 


N  =  PN, 


'0  / 


2A?^?  a2 

e  C  C  d(£-) 

c  c 


(B4) 


which  agrees  with  Equation  (52)  of  Reference  1,  thus  showing  the 
correspondence  of  the  approach  used  herein  with  that  used  in 
Reference  1. 

Equations  (B2)  or  (b4)  imply  then  a  universal  curve  relating 
N/PNq  to  tl/Aoq.  (Henceforth  A  is  to  be  regarded  as  the  compo¬ 
site  value  A  ,  but  the  subscript  will  be  dropped  for  brevity  of 

V 

notation.)  The  consideration  of  much  operational  flight  data  in 
this  light  has  led  to  the  tentative  universal  curve  shown  in 
Figure  1.  It  is  found  that  in  the  range  of  x/Ao  >  10,  the  curve 

V 

is  expressed  well  by 


N 

PN0 


(B5a) 
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and  that  for  x/A a  <  10,  the  relation 

C 


N 

PN0 


b  x 
ae  acc 


(B5b) 


applies  with  reasonable  accuracy.  The  joining  point,  with  equal 
height  and  slope  values  is  given  by 


For  the  curve  shown  in  Figure  1,  values  of  a  and  b  are  found 
to  be  a  =  3.36,  b  =  6.74,  so  that 


(B6a) 


=  e 


x 

Acc 


<  9.13 


(B6b) 


With  A  and  N0  given  for  a  particular  airplane,  and  a  and  P 
stipulated,  the  curve  of  Figure  1  defines  completely  the  expected 
number  of  times  per  second  a  given  response  level  x  is  crossed. 
The  number  of  times  n  the  level  x  is  crossed  in  the  total 
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flight  of  the  aircraft  follows  simply  as 


n  =  W2 

For  N  given  by  Figure  1  and  specifically  by  Equation  (Bfia), 
this  number  becomes 


n 


/3*36 
■  m°(— 7 


(W) 


This  equation,  or  Figure  1,  thus  embodies  a  concise  and  succinct 
statement  for  the  design  of  aircraft  due  to  gust  encounter,  as  is 
discussed  in  the  body  of  the  report.  As  a  reiteration,  the  values 
of  a  =  3.36  and  b  =  6.7-4  are  tentative,  and  should  be  updated 
as  more  and  more  operational  data  are  cast  in  the  '‘’universal” 
form  suggested  here. 


APPENDIX  G 


FATIGUE  CONSIDERATIONS 

Recent  fatigue  studies  v.:'  »■  \ally  loaded  specimens  under 
random  loading  (Reference  3)  indicate  that  the  fatigue  life 
appears  to  be  governed  almost  wholly  by  the  r.m.s.  value  of  the 
applied  stress,  independent  of  the  spectral  shape  of  the  loading. 
Typical  results  are  shown  in  Figure  7.  These  results  are  consi¬ 
dered  significant  and  remarkable  and  suggest  a  new  way  of 
considering  the  fatigue  life  of  aircraft  structures.  Instead  of 
considering  different  loading  cycles  of  different  magnitudes  or 
of  trying  to  considei  fatigue  through  use  of  cumulative  damage 
concepts,  we  advance  here  the  following  hypothesis:  for  rardom 
loading,  as  encountered  by  aircraft,  fatigue  life  Is  a  function 
simply  of  the  r.m.s,  stress  level.  The  r.m.s.  level  may  be 
changing  with  time,  as  large  or  low  level  load  sequences  are 
encountered,  but  it  is  the  over-all  r.m.s.  level  that  is  effective 
and  important.  Witn  a  llucie  stretch  in  thought,  it  is  conceivable 
that  fatigue  due  to  all  loads  sources,  whether  taxiing,  ground-to- 
air  cycle,  maneuvering,  or  gusts  may  be  treated  in  this  way;  hence 
the  r.m.s.  stress  value  of  the  lumped  time  history  of  all  these 
load  sources  is  postulated  as  being  the  primary  value  governing 
fatigue  life. 

The  results  shown  in  Figure  7  are  fitted  well  by  the  expression 

no  -  (I2^0)5'27 

where  n^  is  the  total  number  of  zero  crossings  with  positive 
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slope  necessary  to  produce  fatigue  failure.  The  equation  suggests 
the  general  form 


(Cl) 


where  sr  is  some  empirical  reference  stress.  Consider  that  this 
reference  stress  is  expressed  in  terms  of  the  design  stress  level 
x  according  to  the  relation 


sr  =  fix  (C2) 

and  recall  that  from  Appendix  A  a  is  given  by  (with  a  replaced 

X  W 

by  the  composite  value  c  ) 


Equation  (Cl)  would  thus  appear 


(C3) 


In  application  to  aircraft  flight  through  turbulence,  nQ  would 
be  given  by  nQ  =  NQTt  =  PTNqJ  hence  Equation  (C3)  becomes 


PTN0  = 


(C4) 


The  remarkable  correspondence  between  this  equation  for  fatigue 
and  Equation  (B7)  is  dealt  with  in  the  body  of  the  report. 
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